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Abstract 

Elliptic flow at RHIC is computed event-by-event with NeXSPheRIO. We show that when sym- 
metry of the particle distribution in relation to the reaction plane is assumed, as usually done in 
the experimental extraction of elliptic flow, there is a disagreement between the true and recon- 
structed elliptic flows (15-30% for r?=0, 30% for p±=0.5 GeV). We suggest a possible way to take 
into account the asymmetry and get good agreement between these elliptic flows. 
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Hydrodynamics is one of the main tools to study the collective flow in high-energy nuclear 
collisions. Here we discuss results on elliptic flow obtained with the hydrodynamical code 
NeXSPheRIO. It is a junction of two codes: NeXus and SPheRIO. The SPheRIO code is 
used to compute the hydrodynamical evolution. It is based on Smoothed Particle Hydro- 
dynamics, a method originally developed in astrophysics and adapted to relativistic heavy 
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ion collisions Its main advantage is that any geometry in the initial conditions can be 
incorporated. The NeXus code is used to compute the initial conditions T^, j M and on 
a proper time hypersurface Q. NeXSPheRIO is run many times, corresponding to many 
different events or initial conditions. At the end, an average over final results is performed. 
This mimicks experimental conditions. This is different from the canonical approach in hy- 
drodynamics where initial conditions are adjusted to reproduce some selected data and are 
very smooth. This code has been used to study a range of problems concerning relativistic 
nuclear collisions: effect of fluctuating initial conditions on particle distributions y| , energy 
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dependence of the kaon effective temperature |4|], interferometry at RHIC [5], transverse 
mass distributions at SPS for strange and non-strange particles jfj, effect of the different 
theoretical and experimental binnings , effect of the nature of the quark-hadron transition 
and of the particle emission mechanism jg]. Here we study the evaluation of elliptic flow 
using the so-called standard method. The version of NeXSPheRIO used here has a first 
order quark-hadron transition, sudden freeze out and no strangeness conservation. The only 
parameter, the freeze out temperature, was assumed to be 150 MeV, since this gives good 
agreement for the charged particle pseudo-rapidity and transverse momentum distributions 
for all PHOBOS centrality windows. 

Theoretically, the impact parameter b is known and varies in the range of the centrality 
window chosen. The theoretical, or true, elliptic flow parameter at a given pseudo-rapidity 
r] is defined as 

b J d 2 N/d(j)dr l cos[2((j) - (j) b )}d(t) 

<wJ>-< jdw/dwndcf) > [i) 

(fit, is the angle between b and some fixed reference axis. The average is performed over all 
events in the centrality bin. 

Experimentally, the impact parameter angle <pb is not known. An approximation, ip2, is 
estimated. Elliptic flow parameter with respect to this angle, v^irf), is calculated. Then a 
correction is applied to v^irj) to account for the reaction plane resolution, leading to the 
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experimentally reconstructed elliptic flow parameter v^irj). For example in a Phobos-like 
way 0] 

< v r 2 ec ( V ) >=< <S{V) > (2) 

V< cos[2(^<° - ^> )] > 

where 

obs { \ _ E» (PN/d&drjcos^i - j) 2 )} 

and 

In the hit-based method, and f/^ are determined for subevents rj < and > 
respectively and if v<i is computed for a positive (negative) rj, the sums in 7/> 2 , eq. (4), are 
over particles with rj < (rj > 0). In the track-based method, ip2° an d are determined 
for subevents 2.05 <| r\ |< 3.2, the sums in ip 2 , eq. (4), are over particles in both sub-events, 
t>2 is obtained for particles around < r] < 1.8 and reflected (to account for the different 
multiplicities between a subevent and the sums in eq. (4), there is also an additional y/2a 
with a ~ 1, in the reaction plane correction in eq. (2)). Since both methods are in agreement 
but only the hit-based method covers a large pseudo-rapidity interval, we use this latter 
method. 

We want to check whether the theoretical and experimental estimates are in agreement, 
i.e., < v^irj) >=< V2 ec (r]) >. A necessary condition for this, from eq. (2), is, < v\(rj) >>< 
v% bs {ri) >. In figure 1, we show the results for < v^f]) > (solid line) and < v^irj) > 
(dashed line). We see that < v^T}) ><< v^ s {r}) > for most r/s. So, as shown also in the 
figure, dividing by a cosine to get < f£ ec (?7) > (dotted curve) makes the disagreement worse: 
< 17^(77) > and < v r 2 ec {j]) > are different. This is true for all three Phobos centrality windows 
and more pronounced in the most central window. 



Since the standard way to include the correction for the reaction plane resolution (eq. 
(2)) seems inapplicable, we need to understand why. When we look at the distribution 
d 2 N/d(pdr] obtained in a NeXSPheRIO event (presumably also in a true event), it is not 
symmetric with respect to the reaction plane. (We recall that the reaction plane is the plane 
defined by the impact parameter vector and the beam axis.) This happens because i) the 
incident nuclei have a granular structure, ii) the number of produced particles is finite. The 
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FIG. 1: Comparison between various ways of computing V2 using NeXSPheRIO for Phobos 15-25% 
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centrality window 9]: solid line is v 2 , obtained using the known impact parameter angle 4>b, dashed 
(dotted) line is v 2 bs (v 2 ec ), obtained using the reconstructed impact parameter angle fa without 
(with) reaction plane correction. 

symmetry might be better with respect to the plane with inclination fa in relation to the 
reference axis and containing the beam axis. Therefore we must write for each event 



— - = v fa)[l + £2^(77) cos(n(0-&)) + £ 2< 6 (?7)si n (n(0-0 6 ))] (5) 
acpcLrj n n 

= v (t)) [1 + E 2v^( V ) cos(n(0 - fa)) + £ 2 <° bS (v) sm(n(0 - fa))) (6) 

n n 

It follows that 

v°2 bs (v) = v b 2 (v) cos[2(^ 2 - fa)] + v' 2 b (rj) sm[2{fa - fa)} (7) 

We see that due to the sine term, we can indeed have < v 2 {rj) >>< v\{rj) >, and therefore 
< v r 2 ec ij]) >>< ^2(77) > as in figure 1. The sine term does not vanish upon averaging on 
events because if a choice such as eq. (4) is done for fa, v 2 b (r]) and sin(2(-0 2 — 4>b)) have 
the same sign. This can be visualized with fig. 2a. If the momentum distribution, instead 
of being symmetric with respect to the reaction plane, (for example v\ > 0, v' 2 b = 0) has 
a positive sine term added (v' 2 b > 0), it now points at an angle between and 7r/4 above 
the reaction plane. This angle is in fact ip2 and is determined experimentally with eq. 
(4). Therefore v 2 sin(2(^>2 — <Pb)) > 0. Similarly, if v 2 < 0, ip2 is between — 7r/4 and and 
v 2 sm(2(ip 2 — (f)i ) )) > 0. (Rigorously, this sign condition is true if fa is computed for the same 
77 as v 2 (r]). Due to the actual way of experimentally extracting fa, we expect this condition 
is approximately satisfied only for particles with small or moderate pseudorapidity, which 
are close enough to where fa was computed.) 
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FIG. 2: Assuming (top) d 2 N/d(f)dr] = 1 + 2v\ cos(2(0- <p b )) + 2vf sin(2(0 - <j> b )) with 6 = 0: dash- 
dotted momentum distribution is symmetric with respect to the reaction plane (v 2 > 0, v 2 b = 0) 
and solid is asymmetric [v\ > 0, v' 2 b > 0); assuming (bottom) d 2 N/d<f)di] = 1 + 2v l 2 S cos(2(0 — 
V>2)) + 2v 2 obs sin(2(</> — ^2)) with (fit, = 0: dashed momentum distribution is symmetric with to the 
plane inclination ip2 above the impact parameter and containing the beam axis (v 2 bs > 0, v 2 obs = 0) 
and solid is asymmetric (v 2 bs > 0, v 2 obs > 0). 

In the standard approach, for example as in Phobos analysis, it is assumed that d 2 N/ d(pdr] 
is symmetric with respect to the reaction plane and there are no sine terms in the Fourier 
decomposition in (eq. (5)); eq. (7) leads to (for the hit-based or track-based method) 

< v b (r)) >=< vtiv) > I < cos[2(^ 2 - <j> b )] > (8) 

Then using < cos [2 (ip 2 -0b)] >=< cos[2(ip^ — <j) b )] >=< cos[2(ip^ -<p b )} > and < cos[2(ip% — 
1^2)] >=< cos[2(ip2 ~ 4>b)] >< cos [2 (1^2 — 0b)] >=< cos [2 — 0b)] > 2 (where it is assumed 
that the distributions of ip^ ~ 0b an d i>2 ~ 0b are symmetrical with respect to the reference 
axis and tp^ ~ 0b an d i>2 ~ 0b are independent), eq. (2) follows. However as explained above, 
the use of NeXus initial conditions leads to d 2 N/d4>dr] not symmetric with respect to the 
reaction plane (and presumably this is also the case in each real event), so eq. (8) and (2) 
are not valid. 

As already mentioned, the symmetry might be better with respect to the plane with 
inclination ip 2 i n relation to the reference axis and containing the beam axis. From (5) and 
(6), we have 

v b 2 (v) = vtiv) x cos[2(^ 2 - 6 )] + v' 2 obs ( V ) x sin[2(^ 2 - b )]. (9) 

If the symmetry is perfect v' 2 obs = 0. Otherwise, looking at fig. 2b, if the angular distribution, 
instead of being symmetric with respect to the axis with inclination ip 2 in relation to the 
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impact parameter, (for example v 2 bs > 0, v' 2 = 0) has a positive sine term added (v' 2 > 0), 
it now points at an angle i/j 2 ew greater than ijj 2 . If a negative sine term is added (v' 2 obs < 0), 
it now points at an angle ip 2 ew smaller than ip 2 . Both possibilities are equally likely for a 
given ip 2 but lead to opposite signs for v' 2 obs (r]) x sin[2(-?/>2 ew ~ (fib) ( m general). Therefore 

< v 2 obs (r]) x sin[2(^ 2 — (fib)] >— 0. So whether the symmetry is perfect or approximate, 

< v^irj) >~< v 2 bs (rj) x cos [2 (^2 — </>&)] > and instead of eq. (2) we would have 

< v* ec {ri) >=< v° 2 bs {r]) x y/< cos[2(^ 2 <0 - t/> 2 >0 )] > > (10) 

In figure 3, we show < v 2 ec (rj) > (dash-dotted line) and < v 2 (r]) > (solid line). We see that 
the agreement between both methods is improved compared to figure 1. 



We have also computed the elliptic flow parameter as function of transverse momentum 
for charged hadrons with < rj < 1.5 for the 50% most central collisions. We found that 
< v 2(P-0 > computed as in eq. (1) is well approximated by < V2 ec (p±) > computed as in 
eq. (10). 

In summary, from figure 1, elliptic flow estimated from the standard method with reaction 
plane correction is an overestimate of true elliptic flow (t>2 ec > v 2 ). From figure 3, using a 
method that takes into account the more symmetrical nature of particle distribution in 
relation to the plane with inclination ip 2 with respect to the reference axis and containing 
the beam axis (rather than with respect to the reaction plane), we get a better agreement 



between reconstructed and true elliptic flows (v 2 ~ v 
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0. 



As for overestimating the true elliptic flow, a similar conclusion was reached in [12] and 



In 



121 ] . elliptic flow was assumed proportional to eccentricity and eccentricity was 



computed event-by-event using a Monte Carlo Glauber calculation. As in our case, b is 
known. It was found that the integrated true v 2 is smaller than v r 2 c computed with a two- 
particle cumulant method (for all centralities) and larger than v 2 ec computed with higher 
order cumulants (for centralities 0-80%). In 13j, elliptic flow was computed event-by-event 
within the UrQMD model. Again b is known. It was found that the integrated true v \ is 
smaller than v 2 ec computed with a two-particle cumulant method (for all centralities) and 
equal to v 2 ec computed with higher order cumulants (for centralities 10-50%). Differential 
elliptic flow was also studied leading to similar conclusions. 
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(25 - 50) % 




FIG. 3: Comparison between true elliptic flow v\ (solid line) and suggested method to compute 
reconstructed elliptic flow from data v^ ec (dash-dotted) for the three Phobos centrality windows ^]. 
Squares represent Phobos data (black error bars are 1 a statistical errors and grey bands, systematic 
uncertainties at ~90% confidence level). 



In these two papers, it is expected |l2|,[l3| that there will be differences between v\ and v™ 
calculated with the reaction plane method or two-particle cumulant method both because 
of the so-called non-flow correlations (overall momentum conservation, resonance decays, 
jet production,etc) and event-by-event fluctuations (mostly eccentricity fluctuations). In 
principle, higher-order cumulant methods take care of non-flow effects. If there is still 
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disagreement between the true elliptic flow and higher-order cumulant methods, as in 
then fluctuations are important. If there is agreement as in then non-flow effects are 
important and not fluctuations. In addition to the disagreement between their conclusions, 



12| and do not (neither are expected to) reproduce the RHIC data. So an interesting 
question is whether a more accepted hydrodynamical description would lead to a sizable 
effect. Using NeXSPheRIO, we found that true elliptic flow v\{j] = 0) is overestimated by 
~ 15-30 % (according to centrality) with the reaction plane method, and v\{p±) by ~ 30% 
at p_i_=0.5 GeV. In our case, since < v\ >~< v^ ec >, a large part of the difference between 
the true < v \ > and reconstructed < v^ > is due to the (wrong) assumption of symmetry 
of the particle distribution around the reaction plane, made to obtain < >. 

Finally, we would like to emphasize that it is important to have precise experimental 
determination of elliptic flow, in particular free from the assumption of symmetry that we 
discussed. Elliptic flow teaches us about the initial conditions and thermalization, in princi- 
ple. In this manner, in Q], the author showed that with his hydrodynamical code plus freeze 
out, he could not reproduce 1)2(11), in particular at large 77, and therefore concluded that there 
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might be a lack of thermalization for these large 77's. In [15], it was shown that agreement 
with 1*2(77) data could be obtained for central collisions with a similar hydrodynamical code 
but with color glass initial conditions and, instead of freeze out, a transport code matched 
to the hydrodynamical code to describe particle emission. It was therefore concluded that 
some viscosity was necessary in the hadronic phase. Lastly, in (see figures 3 and 4), it 
was shown that for all centralities, Glauber-type initial conditions plus hadronic dissipation 
lead to a reasonable agreement with data while color glass condensate initial conditions 
plus hadronic dissipation do not, except in the most central window (unless some additional 
dissipation occurs in the early quark gluon plasma phase). Both sets of initial conditions 
without hadronic dissipation tend to underestimate 1*2(77) data if 7/. oni = 169MeV and over- 
predict them if Tf out = lOOMeV. However these conclusions would be affected if 1*2(77) data 
were lower, as we think they should be. (Incidently, though our objective was not to repro- 
duce data, note that our model with freeze out (no transport code) reproduces reasonably 
both the 1*2(77) data as in [lf| (figure 3) and the V2(p±) data (not shown).) Therefore, to 
know e.g. what the initial conditions are or if there is viscosity and in what phase, we need 
to settle the question of whether event-by-event fluctuations are important and take them 
into account in the experimental analysis. 
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